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Abstract. We discuss how adiabatic potentials can be used to create addressable 
lattices on a subwavelength scale, which can be used as a tool for local operations and 
readout within a lattice substructure, while taking advantage of the faster timescales 
and higher energy and temperature scales determined by the shorter lattice spacing. 
For alkaline-earth-like atoms with non-zero nuclear spin, these potentials can be made 
state dependent, for which we give specific examples with 171 Yb atoms. We discuss in 
detail the limitations in generating the lattice potentials, in particular non-adiabatic 
losses, and show that the loss rates can always be made exponentially small by 
increasing the laser power. 
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1. Introduction 

Optical lattices in 3D generated by counterpropagating laser fields provide a periodic 
array of microtraps for cold atoms, and implement quantum lattice models governed 
by Hubbard Hamiltonians flj. This is of interest both in the study of strongly 
correlated systems from condensed matter physics, and in the implementation of 
quantum information processing. Two key elements of cold atom dynamics in 3D 
optical lattices are: (i) The lattice spacing is bounded below by the wavelength a = A/2. 
Therefore, all the relevant time, energy and temperature scales are determined by the 
recoil energy, Er = h 2 k 2 /(2m), where k = 2-k/X is the wavenumber, and m is the the 
mass of the atom. In particular, the hopping matrix element between neighbouring 
sites within the region of validity of the single-band Hubbard model is limited by 
J <C Er. This limits all corresponding timescales, for example the timescale over 
which entanglement can be generated using exchange interactions [21 [3] . (ii) Varying 
the laser parameters provides global control of the lattice potential, but development of 
tools for local operations and readout remains a major challenge in light of applications 
in quantum computing and quantum simulation [4]. Here we show how adiabatic 
potentials can be used to address these issues by creating addressable, state-dependent 
sub-wavelength lattices. Each lattice period can be divided into a set of wells, which 
can be frequency-addressed via their offset energies, reminiscent of quantum computing 
schemes with electrons in arrays of electrically gated quantum dots [5]. 

In recent seminal experiments [21 E], optical superlattices have been used to provide 
addressability and to entangle atoms in a series of double-wells in one parallel operation. 
Our goal here is to achieve similar control on sub-wavelength scales by coupling optical 
lattice potentials with dressing fields, with the benefit that all corresponding timescales 
are significantly faster, and energy and temperature scales significantly higher. Single 
particle effects of such dressed potentials j6j [7J El El Ell EU EI21 EH [H] have been studied, 
e.g., in the context of lithography [15], laser cooling [161 E]) an d Radio Frequency 
(RF)-dressed potentials on atom chips [18]. Here we show how the adiabatic potentials, 
created by coupling optical lattice potentials with additional dressing fields, can be used 
for subwavelength control in the context of manybody physics. We will demonstrate how 
adiabatic potentials can be used to generate state-dependent, addressable subwavelength 
lattices. In such a subwavelength lattice, each lattice period is subdivided into several 
potential wells, for which the relative well depth can be tuned by either applying external 
magnetic or electric fields, or by adjusting the laser parameters, e.g. the Rabi-frequency 
and detuning. Atoms loaded into different wells within the subwavelength structure can 
also be frequency-addressed via their different offset energies. 

In view of recent seminal progress with alkaline-earth-like atoms, e.g. in the context 
of developing optical clocks [191 EQ]> in many-body cold atom physics [21], and in 
applications in quantum information [221123], we investigate, in particular, the possibility 
for generating sub-wavelength potentials with these species. The existence of metastable 
triplet states makes possible the production of near-resonant optical lattices, which allow 
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Figure 1. Sub- wavelength lattice with off-resonant and near-resonant A/2 background 
potentials (see text), (a) Schematic of a controllable subwavelength lattice, where 
each lattice period of an initial A/2 is subdivided into several potential wells whose 
position and well depth can be controlled by changing the laser parameters. Note that 
addressability in the subwavelength lattice is modulo A/2, with A the wavelength of 
the laser generating the background lattice potentials, (b) (left) Level structure and 
laser configurations for generating background potentials V g . e (x) in the off- resonant 
scheme, (right) Dressing fields with Rabi-frequency — 1,2,3,...) and detuning 
Si are applied to couple the background potentials to generate subwavelength lattice 
V^{x). The separation between neighbouring wells is approximately A/2(n + 1) with n 
dressing fields, (c) Level structure and laser configurations for generating near-resonant 
subwavelength potentials in a two level system. The counteroscillating background 
potentials V^(x) are now induced by directly coupling the two levels \g) and |e) with 
standing wave fields with Rabi-frequency flo(x) and detuning (5o(red). Dressing fields 
with Rabi-frequency OA* = 1,2, ...) couple the background potentials to generate the 
subwavelength potential. 



for the creation of state-dependent sub-wavelength lattices. We show how this can be 
used to generate and manipulate entanglement, making use of the spin dependence 
and/or shorter tunneling times within the subwavelength structure. Of course, in the 
context of adiabatic potentials, the question of non-adiabatic losses will always arise. 
We analyse these losses in detail in the present context, and conclude that they can 
always be suppressed at the expense of using additional laser power. 

This paper is organised as follows: We first give a general description in Sec. [2] of 
two schemes to create adiabatic potentials with two-level atoms. In Sec. [3j we extend the 
idea to multi-level atoms, describing the additional control that is present when these 
lattices are generated for real atoms. In particular we discuss possibilities for state- 
dependent potentials with alkaline-earth-like atoms. In Sec. HI we add interactions 
between atoms and give examples for possible applications of subwavelength lattices 
in quantum simulation and quantum information. This includes a discussion of how 
entanglement within the subwavelength structure can be generated and controlled. In 
Sec. |5l we analyse the Landau-Zener-type loss rate from the adiabatic potentials in detail 
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and discuss how it should be suppressed. In appendix A we give additional details on 
the clock transition in alkaline-earth-like atoms, which are essential for the generation 
of subwavelength lattices with these atoms. 

2. Generating subwavelength lattice potentials 

In this section we introduce the basic ideas behind addressable subwavelength potentials 
for two-level atoms. This will be generalised to multi-level atoms in Sec. EJ where we 
also discuss their implementation with alkaline-earth-like and alkali atoms, and how 
they can be made internal state- (or qubit-) dependent. We discuss these ideas on the 
single-particle level initially, with interactions between atoms to be added in Sec. HI 

As shown in Fig. la, our goal is to divide each period of a standard A/2 optical 
lattice into several potential wells with controllable barrier heights and well depths in 
this subwavelength structure. A laser can then be used to address atoms on different 
sites independently based on energy selectivity, coupling to an unperturbed reference 
state. In addition, we can control in this way the tunneling of atoms between adjacent 
wells in the subwavelength structure. 

The simplest way to achieve this goal is shown in Fig. lb. We consider a two-level 
atom with two long-lived internal states |e) and \g). We assume that we can generate 
counter-oscillating off-resonant potentials for these two internal states, for which the 
implemenation will be discussed in Sec. [3j We then couple the internal states with 
additional dressing fields, which generate avoided crossings in the background potentials, 
and thus a subwavelength structure. Here, the well depth and barrier heights can be 
controlled by varying the detunings and Rabi frequencies of the dressing fields, which 
shifts the locations of the avoided crossings. In Fig. lb, we show the simplest case with 
isolated avoided crossings. 

It is also possible to generate background potentials with near-resonant coupling 
between \g) and |e) with a Rabi frequency varying sinusoidally in space (e.g., as 
generated by a standing wave light field). This results in counter-oscillating potentials 
due to AC-Stark splitting, which can then be coupled by a series of dressing fields to 
produce the desired subwavelength structure (see Fig. lc). This scheme is motivated, 
in particular, by alkaline-earth-like atoms (as discussed in Sec. [3]). 

2.1. Addressable subwavelength lattices with two-level atoms in ID: Off-resonant 



We now discuss the details of the generation of subwavelength lattices with off-resonant 
background potentials for a two-level atom. In the rotating wave approximation with 
respect to the optical frequencies, we may write the single atom Hamiltonian as: 



scheme 



H = H M + H + Hi, 



(1) 



with 




(2) 
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#o = V g (x)\g){g\ + (V e (x)-5 1 )\e)(e\ 

Hi = Yl ^exp(-t(n-l)5t)\e)(g\+h.c.. 

n=l,2,... 

Here, H M is the kinetic energy, H describes the background potentials, and H 1 contains 
the dressing fields with Rabi frequencies Q n . For dressing fields produced as an equally- 
spaced comb of sidebands, the detunings can be written as S n = Si + (n — 1)5, where 
Si is the detuning of the first dressing field. In the off-resonant scheme, the background 
potentials for the two internal states are generated by standing wave fields and are 
opposite to one another. For this, we may take V g (x) = — V e (x) = —V sm 2 (kx), where 
k = 2n/\ is the wavenumber of the laser generating the background potential. 

Generation of the adiabatic potentials is based on the validity of a Born- 
Oppenheimer-type assumption. This involves the fact that the kinetic energy of the 
atoms is small on a scale given by the separation of the resulting adiabatic potentials. 
The wavefunction \&(x, t)) of a single atom satisfies the Schrodinger equation, 

ih^Mx, t)) = (H M + H + Hi)mx, t)). (3) 

If we omit the kinetic energy term from the Hamiltonian, we obtain an equation 
for adiabatic eigenstates, \^ a (x,t)), 

ih ^_ e -iE a (x)t^ a) = {H(j + Hl ) e -iE a (x)t^ a ^ (4) 

where the \^ a (x, t)) are periodic functions and can be expanded as a Fourier series. For 
a two-level atom, it takes the form 

oo 

\* a (x,t)) = exp(-USt)[c l g (x)\g) + c l e (x)\e)]. (5) 

l=— oo 

This results in a Floquet eigenvalue equation for E a (x), which plays the role of the 
adiabatic potentials, 

(V g (x) - lS)c l g + \^ n - X = E a (x)c l g (6) 

n 

(V e (x) -Si- lS)c{ + J2 Y° l 9~ n+1 = E ^ C - W 

n 

The complete wavefunction can then be expanded in a basis of these adiabatic 
eigenstates, which play the role of Born-Oppenheimer channel functions, 

\<f>(x,t)) = Y,c a (x,t)\V a {x,t)), (8) 

a 

resulting in the equation 
d_ 



ih—c a (x, t) = [H M + E a (x)]c a (x, t) + Y H M c pi x , t), (9) 



P 

a/3 



where = (i&a\HM\^p)- The equations for different a decouple provided that the 
mixing terms coupling the channels are small. We expect this to be the case provided 
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Figure 2. Illustration of the controllability of subwavelength potentials, (a) Typical 
subwavelength lattice with n = 1. The potential in the figure is generated with 
the off-resonant scheme with the parameters: Vq = 200kHz, ill = 2ir x 50kHz, 
S\ = 120kHz. Similar potentials can also be generated with the near-resonant scheme, 
(b) Subwavelength double-well structure generated with the off-resonant scheme with 
the parameters: V = 200kHz, Qi = 2tt x 50kHz, ft = 120kHz, 2 = 2vr x 50kHz, 
S = 220kHz. The relative depth of the wells in the subwavelength structure can be 
tuned by changing f?2 or 5. (c) Subwavelength double-well structure generated with 
the near-resonant scheme with the parameters: Oq = 27r x 400kHz, So = 100kHz, 
S7i = 2tt x 560kHz, S — —400kHz. The central barrier height in the lower dressed 
potential can be adjusted by changing either S or f2i. (d) Subwavelength double- 
well structure generated with the off-resonant scheme with n = 3 dressing fields: 
Vq = 200kHz, Si = 60kHz, fli = fl 2 = fl 3 = 2tt x 50kHz, 5 = 140kHz. 



that the kinetic energy of the atom in the lattice, which is typically on the order of 
the recoil energy, is much smaller than the separation between the adiabatic potentials, 
which is characterised by laser parameters, e.g. the Rabi frequency and detuning. Below 
we will initially neglect these terms, and in Sec. [5] we will discuss in detail non-adiabatic 
processes, and how they should be suppressed in experimental implementations. 

In the case of two-level atoms, the eigenvalues E a (x) of the equations appear in 
pairs, and may be identified as the upper V + (x) and the lower V~(x) dressed potentials 
for a given Floquet manifold (see Fig. 2). The corresponding eigenstates for the adiabatic 
potentials in general have components from different Floquet manifolds. For the simple 
case of one dressing field in the off- resonant scheme (with only Q\ non-zero), the different 
manifolds in the Floquet basis decouple, and we only need to diagonalise a two-by-two 
matrix. In this case, the expression for the adiabatic potential and its corresponding 
states are: 

V ± = - - ± VX 2 + Y 2 

|*±) = I[(l T - X )h g ) ± (l ± - X )h e )V (10) 
1 ' 2 U y/X 2 + Y 2 VX 2 + Y 2 

where X = V e (x) — 5i/2, Y = Q\/2, Z = 5\. The resultant dressed potential has a lattice 

spacing of A/4. The structure of the adiabatic potential becomes more complicated with 
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more dressing fields. With n dressing fields, it is possible to create a lattice with spacing 
approximately A/2(n + 1). This increase in the spatial resolution is at the expense of a 
decrease in the lattice depth. 

2.2. Addressable subwavelength lattices with two-level atoms in ID: Near-resonant 
scheme 

In the near-resonant scheme, the counter-oscillating background potentials are generated 
by directly coupling the two internal states with standing wave fields, producing a pair 
of AC-Stark split states. The Hamiltonian describing the coupling between states is 
now given by [cf. Eq. (j2J)] 

H = -5o|e><e|+(^|e><s|+h.c.) 

Hl= J2^eM-™6t)\e)(g\ +h.c, (11) 

n 

where flo(x) = Q sin(A;x) is the Rabi frequency corresponding to a standing wave field, 
and Q n (n = 1,2,...) are the Rabi frequencies of the dressing fields with corresponding 
detunings 5 n = 5q + n5. Note that the background potentials are associated with the 
eigenstates of H , with the values of the the potentials and the states given by the same 
form as Eq. §M$, with X = -5 /2, Y = Vt (x)/2 and Z = 5 [Uj. 

Similar to the procedure in the previous section, we derive the equations for the 
adiabatic potentials, 

n 

(So - m + + y: Y cl r = E °w<t» ^ 

n 

with E a (x) being the adiabatic potentials. 

In Fig. 2 we give examples of the flexibility in constructing adiabatic potentials via 
the near-resonant scheme as well as the off-resonant scheme, including the realisation 
of double-well structures on a subwavelength scale. In Fig. 2b, c, we demonstrate the 
simple case of one dressing field: by varying the Rabi-frequency or the detuning of the 
dressing field, the spatial structure and the central barrier height of the subwavelength 
double-well structure can be controlled. This is essential for entanglement manipulation 
and state preparation, as we will discuss later. 

2.3. Adiabatic potentials in higher dimensions 

In higher dimensions, adiabatic potentials can also be used both to create shorter 
period lattices, and to change the geometry of the lattices. Two examples of 2D 
potentials are shown in Figs. El HI In Fig. [3] we show the form of the dressed potentials 
obtained with an initial square lattice with period a. From initial off-resonant potentials 
V e (x,y) = Vosin 2 (A;;x) + Vq sin 2 (kiy) , and V g (x, y) = —V e (x,y), and using a detuning 
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Figure 3. Adiabatic dressed potentials in 2D from an initial square lattice, with depth 
2Vo, formed with a single dressing field, Si — 2V (Note that the potential energy 
represented by each colour is particular to the figure part, and in each case darker 
colours correspond to lower potential energy) : (a) The initial off-resonant potential for 
state \g), V g (x, y) — — Vq sin 2 (fc;a;) — Vq sin 2 (fc;?/). (b) The initial resonant potential for 
state |e), V e (x, y) = — V g (x, y). (c) The lower dressed potential V~ (x, y), corresponding 
to a square lattice with l/\/2 times the original period, (d) The upper dressed potential 
V + (x,y), corresponding to interconnected rings. Note that the adiabatic potentials 
V ± (x, y) can be found from Eq. (10), and have a maximum depth Vq, which decreases 
as the Rabi frequency of the dressing field increases. 



Si = 2Vo, we obtain for V~(x, y) a square lattice of period a/\/2, as the minima in this 
dressed potential occur at each position of a minimum in V e (x, y) or V g (x, y). In the 
upper adiabatic potential V + (x), we observe the same form inverted, giving rise to a 
potential consisting of interconnected rings. 

In Fig. HI we show the dressed potentials arising from an initial triangular lattice, 
formed by three interfering laser beams at an angle of ir/3 to each other. This produces 
a well-defined Honeycomb lattice structure, with auxiliary sites introduced in the lower 
dressed potential, V~(x, y). In the upper dressed potential, V + (x,y) we observe a 
series of disconnected rings (see Fig. HH). Such disconnected ring structures can also be 
produced from a square 2D lattice with two dressing frequencies. Note that the control 
over well depth here could have many useful applications, e.g., in the case of the square 
lattice it could be well-suited to implementation of ring-exchange Hamiltonians [21] . 



3. State-dependent subwavelength potentials with multi-level atoms 

The ideas for creating addressable, subwavelength lattices for two level atoms as 
discussed in the previous section will now be extended to multi-level atoms, where 
it is possible to create internal state- (qubit-) dependent potentials. We will discuss, in 
particular, the case of alkaline-earth-like atoms with nonzero nuclear spin, e.g. 171 Yb. 
In addition to the optical excitation to the metastable triplet levels, we have an extra 
degree of freedom in the form of different magnetic sublevels, which arise from the 
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Figure 4. Adiabatic dressed potentials in 2D from an initial triangular lattice, with 
depth 2Vb, formed with a single dressing field, S\ = 2Vb (Note that the potential energy 
represented by each colour is particular to the figure part, and in each case darker 
colours correspond to lower potential energy) : (a) The initial off-resonant potential for 
state \g), V g (x,y). (b) The initial resonant potential for state |e), V e (x,y) = —V g (x,y). 
(c) The lower dressed potential V~(x,y), corresponding to a honeycomb lattice with 
additional sites in the centre of a ring-shaped barrier, (d) The upper dressed potential 
V + (x,y), corresponding to disconnected rings. Note that the adiabatic potentials 
V^{x, y) can be found from Eq. (10), and have a maximum depth Vo, which decreases 
as the Rabi frequency of the dressing field increases. 



nuclear spin. 

3.1. Off-Resonant Lattices 

3.1.1. Lattice with Alkaline- earth-like atoms As an example, we consider the level 
structure of 171 Yb, for which the nuclear spin / = 1/2. We will make use of the states 
on the clock transition. The states in the ground state manifold are \ l Sa,mi = ±1/2), 
and those in the excited state manifold are | 3 Po, m F = ±1/2). It is also possible to 
generalise the scheme to other alkaline-earth-like atoms with non-zero nuclear spin, e.g. 
87 Sr (/ = 9/2), though one then needs to carefully choose the working states in the 
hyperfine structure [25], [26] . 

To generate off-resonant subwavelength lattices, we induce opposite background 
potentials for the states on the clock transition. This is done by applying a set 
of standing wave laser fields at a specific "anti-magic" wavelength at which the AC 
polarisability of the manifolds 1 So and 3 Po are exactly opposite. This is exactly the 
opposite that used in optical lattice clock experiments [HI EHl E7], where a 

"magic" wavelength is chosen so that the AC Stark shifts of the clock states are equal. 
(We expect that the "anti-magic" wavelength in our case to be ~ 600nm for m Yb, and 
the resultant AC stark shift on the order of several hundred kHz for a laser intensity of 
/ ~ 10kW/cm 2 .) Once the background potentials are generated, we can couple these 
potentials by applying additional optical dressing fields tuned near the resonance of the 
clock transition 1 Sq — 3 Pq at ~ 578nm. Atoms can be loaded into the minima of the 



State- dependent, addressable subwavelength lattices with cold atoms 10 




Figure 5. Atomic level structure and laser configurations for the generation of sub- 
wavelength lattices in Yb from off-resonant background potentials, which are formed 
by AC-Stark shifts represented by the red arrows, (a) Formation of subwavelength 
lattices in two two- level systems with 7r-polarised dressing fields (black arrows), where 
the magnetic field induces differential Zeeman shifts in the 1 S'o (/ii ~ 760Hz/G), and 
the 3 Po (fj.2 ~ 350Hz/G) manifolds 28\. (b) Formation of subwavelength lattices in 
two two-level systems with circularly-polarised dressing fields (black arrows) . 

upper(lower) adiabatic potential by first preparing them in one of the manifolds 1 So or 
3 Po, and then appropriately ramping up the dressing fields adiabatically. 

As shown in Fig. 5a, the and 3 Po manifolds for 171 Yb are each split into two 
sublevels. This gives us a four-level system, and dressing fields coupling the different 
manifolds can be produced using either 7r-polarised or circularly polarised light. In 
appendix A we summarise the origins of the corresponding matrix elements [29, 30J. 
Note, in particular, that the matrix elements for circularly polarised couplings are non- 
zero as a result of contributions in the 3 Po manifold from other levels due to hyperfine 
mixing [29j [30] . In the case that we use only circularly polarised dressing fields, we will 
produce two independent two-level systems, where the resulting subwavelength lattice 
can be manipulated for each by altering the appropriate component of each circular 
polarization in the dressing field. For 7r-polarised light, we also obtain two independent 
two-level systems, in which we can offset the relative well depths by applying an external 
magnetic field. Due to the differential Zeeman shift in the 1 So and 3 Po manifolds 
[28], this will alter the detunings of the dressing fields, as well as the energies of the 
uncoupled internal states. For example, with the laser parameters: Vq = 200kHz, 
Qi = 2tt x 200kHz, 5\ = 200kHz, and under a small magnetic field of B = 10G, the 
estimated energy difference between the potentials corresponding to the two different 
dressed states at kx = 7r/2 is ~ 7kHz, and at kx = tt is ~ 4kHz. These numbers 
increase significantly as the external magnetic field is increased. This allows the state 
to be rotated between dressed states of the different two-level systems selectively for 
different wells of the substructure using RF coupling fields. This would be useful for 



State- dependent, addressable subwavelength lattices with cold atoms 11 

(al) (a2) 




(bl) 

X/2 




Figure 6. Schematics for the generation of an off-resonant subwavelength lattice 
with 87 Rb. (al) Lights of different circular polarisation are tuned near the centre of 
the fine structure splitting to generate spin-dependent potentials. (a2)(bl) RF fields 
are applied to couple different hyperfine states with counter-oscillating potentials, 
which are generated by shifting the spin-dependent potentials. (b2) The resultant 
adiabatic potential in the case of one dressing field has a lattice spacing A/2, with A 
the wavelength of the laser generating the spin-dependent potentials in (al). 

applications to quantum information processing (see below). 

3.1.2. Lattice with Alkali atoms One may also generate off- resonant subwavelength 
lattices with alkali atoms. Here, the counter-oscillating background potentials should 
be produced for two different hyperfine ground states making use of well known ideas 
to produce spin-dependent lattice potentials [311 [32], as shown in Fig. 6. These spin- 
dependent lattice potentials can then be coupled using RF dressing fields to produce 
sub-wavelength optical lattices. Details of this scheme are summarised in the caption 
of Fig. 6. Note that with alkali atoms it is not straightforward to find combinations of 
states that act as independently coupled two-level systems and that are simultaneously 
collisionally stable. This will limit the application of alkali atoms in generating state- 
dependent subwavenlength potentials, as well as in entanglement operations, as we 
will discuss later for alkaline-earth-like atoms. In addition, with alkali atoms, the two 
hyperfine ground states used to generate the dressed potentials decohere fast in the 
presence of magnetic field fluctuations, which is not the case with alkaline-earth-like 
atoms, where the dressed potentials are generated almost purely from nuclear spins. 
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Figure 7. (a) Level structure and laser configurations for the generation of near- 
resonant subwavelength potentials with 171 Yb using cross-coupling dressing fields, (b) 
(above) With only the coupling laser fio, there are two independent sets of spin- 
dependent lattice potentials (i — 1,2), where i labels the two different sets of 
potentials created by lasers with different polarisation and ± labels the upper/lower 
potential in a given set. The relative position of the two sets of potentials can be 
adjusted by varying 9 in the lin-Z-lin configuration. The lattice spacing in this case is 
still A/2, (below) By adding one dressing field Qi, the lower potentials become state- 
dependent subwavelength potentials, with lattice spacing A/4. The relative positions 
of the two lower potentials can be adjusted by varying 9. 

3.2. Near-Resonant Lattices 

3.2.1. Lattice with Alkaline- earth-like atoms The long-lived 3 Po level motivates the 
application of near-resonant lattices, as discussed previously in Sec. I2.2[ with direct 
coupling with a standing wave field tuned to the clock transition producing counter- 
oscillating dressed potentials. The great advantage here, though, is that we can now 
produce these background potentials with lights of two different circular polarisations, 
generating two sets of independent background potentials (see Fig. 7a). We can take 
advantage of this in two ways: (1) We can couple these potentials with circularly 
polarised dressing fields (see Fig. 7), so that we produce two independent sets of 
addressable sub-wavelength lattices; or (2) We can couple these potentials with n- 
polarised fields, so that we form dressed states with all four levels (see Fig. 8a), 
and use this to engineer particular forms for the potentials. For the convenience of 
discussion, we define the notation for the internal states \gi) = | 1 ^o, rrz/ = +1/2), 
\g 2 ) = rSo.mj = - V2>, |ei) = \ 3 Po,m F = -1/2), and \e 2 ) = \ 3 P ,m F = +1/2). 

In the first case, where the states are cross-coupled as shown in Fig. 7, we then 
produce state-dependent lattices in the sense that we have two independent two- 
level systems, in which we form dressed states of \gi) and |ej)(z = 1,2) for each 
i independently. Both the near-resonant background potentials and the final sub- 
wavelength potentials are internal-state dependent, and can be easily manipulated 
independently. For example, using lasers in a linZlin configuration (two counter- 
propagating, linearly polarized beams with a relative polarization angle 8), the lattices 
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Figure 8. (a) Level structure and laser configurations for the generation of near- 
resonant subwavelength potentials with 171 Yb using 7r-polarised dressing fields, (b) 
Atoms in the subwavelength double-well structure can be entangled by a controlled 
evolution with the exchange interaction turned on. The interaction can be switched 
on and off by varying the height of the central barrier, (c) In the presence of a 
small magnetic field and as the angle in the lin-Z-lin configuration 8^0, atoms with 
different internal states in the double-well structure, which can be defined as qubit 
states |0i), (i — L,R) as in the figure, feel different potentials. Atoms within the 
subwavelength structure can thus be selectively addressed. 



can be relatively shifted by varying 9 [3U [32] (see Fig. 7b). Such spin dependence 
could be used similarly to the case in Alkali atoms [3U [32] , and we expect losses due 
to spontaneous emissions to be strongly suppressed in the present case. Note that the 
dipole matrix element for the cross coupling is given by 2ireohc 3 T/(ujQ), where co>o is the 
frequency of the clock transition 1 So — 3 Pq, and r ~ 2tt x lOmHz is the natural linewidth 
of the 3 P state in 171 Yb [271 EHl EHl [30] (see appendix A for more details). 

Subwavelength addressability can be achieved in this configuration in a manner 
similar to that discussed in the previous section. By preparing the atoms near the 
minima of the lower adiabatic potentials and applying an external magnetic field to lift 
the degeneracy of the states in the subwavelength structure, we can perform rotations 
between the independent two-level systems selectively by well within the subwavelength 
structure. 

In the second case we consider using 7r-polarised dressing fields as shown in Fig. 8. 
Here, all the states are coupled to make the adiabatic potentials, and changing the 
angle 9 in the lin-Z-lin configuration modifies the shape of all four adiabatic potentials 
in each Floquet manifold. In this configuration, it is possible to produce a repeated 
double-well potential, in which the two wells can be individually addressed. This is done 
by first setting the angle 9 in the lin-Z-lin configuration to a small but nonzero value. 
This will change the shape of the potentials and lift the symmetry in the double-well 
structure. An example is given in Fig. 9, where at Qq ~ 2n x 400kHz, Q\ ~ 2tt x 480kHz, 

5 100kHz, 5 420kHz, 9 = 0.1, B = 10G, the energy shifts are ~ 1.3kHz and 

~ 11kHz for the two sub- wells, respectively. Atoms in the left and right hand side of 
the double-well structure can then be rotated selectively by a RF coupling after the 




Figure 9. Biased double-well structure generated with the configuration in Fig. 8 
for subwavelength addressing with 171 Yb. The two originally degenerate potentials 
(the dashed and the solid curve) are shifted relative to each other under a small 
bias magnetic field, and with a finite angle 9 in the lin-Z-lin configuration as shown 
in Fig. 8a. Parameters: n Q ~ 2tt x 400kHz, S ~ -100kHz, S ~ -420kHz, 
fix = 2tt x 480kHz, 6 = 0.1, B = 10G. The energy shifts are AE a ~ 1.3kHz, and 
AEf, ~ 11kHz, respectively. 



application of an external magnetic field to shift the potentials. 

In practice, strong coupling at intensity / on the clock transition will also give rise 
to AC-Stark shifts AE of the states | e^) and \gi), which, as AE oc / and f2 oc -\/7, 
will become important at large fields (e.g, ~ 10kW/cm 2 for 171 Yb [33]). These have 
been neglected in the two-level atom picture presented in Sec. 12. 21 At higher fields, 
we essentially obtain a combination of the off-resonant and near-resonant background 
potentials. If the polarisability of 3 Po and 1 So are opposite, it will simply enhance 
the depth of the resulting potentials. These combined adiabatic potentials are still 
spin-dependent as the AC-Stark shifts lose their spatial dependence as the angle 9 is 
increased to 9 = tt/2. In the case that these polarisabilities have the same sign, the 
AC Stark shift can also be cancelled or modified by an additional laser at a different 
wavelength. 

3.2.2. Lattice with Alkali atoms With alkali atoms one could also form near-resonant 
lattices by coupling different hyperfine ground states in the Si level, e.g., via a 
Raman process. If the Raman process beams have a sinusoidally varying effective Rabi 
frequency, this will then generate two counter-oscillating dressed states. These can 
be coupled by RF dressing fields to form subwavelength lattice potentials. However, 
this does not give rise to simple internal state-dependent lattices as in the case with 
alkaline-earth-like atoms. 



State- dependent, addressable subwavelength lattices with cold atoms 



15 



4. Application of subwavelength lattices in quantum simulation and 
quantum information 

Having established the methods to realise subwavelength lattice potentials, we now 
turn to discuss ways to utilise the flexibility in control of the potentials for applications 
in quantum simulation and quantum information. This requires consideration of the 
many-body physics in these potentials. We will show how a generalised Hubbard model 
can be realised by engineering the subwavelength potential. This generalised single- 
band Hubbard model features fast tunnelling rates due to the small spatial separation 
of different wells in the potential, which can be utilised to improve the timescales 
associated with realisation of interesting quantum phases and entanglement via exchange 
interactions. 

4-1- Single-band Hubbard model for the subwavelength lattice 

We begin by considering the simplest case of the off-resonant scheme with a sinusoidal 
lattice in ID, beginning with two background counter-oscillating potentials V e (x) = 
— V g (x) = V sm 2 (kx), as discussed in Sec. 12.11 We then choose the detuning and Rabi 
frequency of n dressing transitions so that the resulting bound states in each well of the 
adiabatic potential have approximately the same energy. In this way a subwavelength 
lattice will be formed in which the lattice spacing scales approximately as X/[2(n + 1)], 
while the depth V n of the potential can be as large as V n ~ Vo/(n + l) and decreases with 
increasing Rabi frequency Q n . Motion can be considered in the tight-binding limit of the 
resulting lattice model provided V n /[{n + l) 2 ] > 5Er, and the dynamics are described 
by the single-band Hubbard model 



Here, c a is the annihilation operator for the bosonic or fermionic atoms, and a specifying 
the lattice site and internal atomic (spin) state, J a p denotes the tunneling rates. We 
denote by A Q the energy offset for the internal state and lattice site specified by a, 
and U a/ 3 are the collisional interaction energies, which follow from our knowledge of the 
potentials and scattering properties. A single-band Hubbard model is valid provided all 
of these energies are smaller than the separation energy to higher bands, which, in the 
case that the lowest bound states in all wells have the same energy, is bounded above 
by the oscillation frequency in a typical well, to = [AdiVoEji/(d1 + Of) 1 / 2 ] 1 / 2 //?,. 

If we choose the detunings of the n dressing transitions so that A a is constant, 
then in ID, J a p = J n 5 a+ i^ can be estimated via a simple rescaling of the recoil energy 
E R = h 2 k 2 /2m, as 



Here, the tunneling rates increase exponentially with increasing n. Note, however, that 
these rates can vary within the substructure, because the wells with n > 1 dressing 




(14) 




(15) 
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Figure 10. Tunneling rate J n in the presence of n dressing fields as a function 
of Vq/Er. The squares are numerically-obtained values for the case n = 0(lowest 
curve) and n = 1 (upper two curves), respectively. The dotted lines are plots of the 
n-dependent scaling function (see Eq. (15) in the text). The lowest curve is for the 
case n = 0, and the higher two curves are for the n = 1 case. 



frequencies will not be exactly symmetric, and this can give rise to small splittings 
within the lowest Bloch band. In Fig. 10 we plot values of J n for varying depths of 
the background potential, and number of dressing transitions n. We will typically have 
onsite interactions only, which can be controlled via the shape of each well and the 
scattering lengths for atoms in the different internal states \e) and \g). This can give 
rise to generalised Hubbard models, where the onsite interaction changes values between 
neighbouring sites. Such models can support non-trivial phases, such as coexisting 
Mott insulating and superfluid phases [34]. Note that when there is only one dressing 
frequency n = 1, the value of U a p is approximately the same as the corresponding onsite 
interaction in the original background lattice potentials, as the curvature of the wells is 
essentially unchanged. 

4-2. Applications to Quantum Information 

As discussed in Sec. 13.21 we have with 171 Yb the possibility to form different addressable 
subwavelength lattices from four internal states. We have different possibilities then to 
encode qubits and generate entanglement depending on our choice of configurations for 
the dressing transitions. 

In the case that we choose circularly polarised dressing transitions (see Fig. 7), we 
obtain two independent two-level systems, for which we can control the subwavelength 
lattices independently (see Sec. 13 . 2j) . We can then encode qubits on the dressed states of 
the two two-level systems. These can be addressed individually within the substructure 
of the subwavelength lattice (i.e., modulo the period of the background potentials, where 
other recent ideas could be applied to provide addressing on the larger length scale [36]). 
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We can then make use of the state-dependent lattice potentials to perform entangling 
operations, e.g., using controlled collisions [3T| [32]. 

Another possibility arises in the case that we use 7r-polarised dressing fields and 
form adiabatic eigenstates combining all four internal states (See Sec. 13.21 and Fig. 8). 
Here it is possible to create entanglement for fermionic atoms such as 171 Yb within 
the double- well structures using exchange interactions [H [37] , where the fast tunneling 
rates J a p between neighbouring wells give rise to faster gate times than are possible 
without subwavelength lattices. This leads to quantum gate schemes in the spirit 
of the Loss-DiVincenzo proposal for quantum dots |40j, which is based on encoding 
qubits on electron spins in electrically gated quantum dots, and performing entangling 
operations using exchange interactions between neighbouring dots by ramping down the 
corresponding potentials. The same techniques can be applied here. We consider the 
two dressed states of our four initial internal levels that have the same dressed potentials, 
and define these to be our qubit states. We then consider the dynamics generated by the 
Hamiltonian in Eq. (114j) . but now with the tunnelling parameters J a p determined by our 
choice of subwavelength lattices with repeated double-well potentials. The tunnelling 
between these wells can be made negligible on the timescale of the experiment, so that 
the dynamics within each double well is independent. We then have a situation very 
similar to that discussed in the analysis of spin-exchange gates, e.g., in Ref. [3], leading 
from the two-site double-well potential to exchange interactions between the spin states 
in each well, given by H SE = Jxd^.a^ where J\ = J 2 (1/Ul + 1/Ur). Here, o~l and o~r 
represent the spin operators in the left and right well respectively, where the spin in 
each well corresponds to two-state system consisting of the two different dressed states 
we are condsidering. The coefficient J\ is determined from the tunneling amplitude 
between sites, J = J a p [see Eq. (fT4l) ]; and the onsite interaction strengths Ul and Ur, 
between two atoms, one in each state, in the left well and right well respectively. 

In the context of electrically gated quantum dots, one of the most challenging 
aspects of implementing quantum computing is to achieve local addressing of the qubits. 
As described above, this is made possible here by the addressability of the individual 
wells within the subwavelength structure, as discussed in Sec. 13.21 There we specifically 
discussed rotation between different qubit states. However, it is also possible to perform 
readout state-selectively by creating an energy offset between the two qubit states by 
introducing an external magnetic field, and then coupling them state-selectively to the 
3 P2 level. This is possible because the 3 P2 level has a very small natural linewidth 
r ~ 27r x 15mHz [29J, so that coupling to this level can be easily made frequency- 
selective. Coupling between the qubit states and 3 P2 can be generated with a Raman 
pulse, making use of the larger coupling elements of both 3 P 2 and 3 Po to high-lying 
states such as 3 Si. 
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5. Landau-Zener-type Losses 



As discussed in Sec. 12.11 one of the major limitations of the adiabatic potentials is 
the loss of atoms via Landau-Zener-type processes. These losses correspond to the 
tunneling of atoms from a given adiabatic potential into the continua of other adiabatic 
potentials. In the following we show that these processes can be made exponentially 
small by increasing the laser power. 

Landau-Zener-type losses can be understood by looking at the simple scenario of a 
single avoided crossing of two adiabatic potentials. This problem has been discussed in 
detail in Ref. [IE] , where the two adiabatic potentials ±U (x) are generated by coupling a 
two-level atomic system using a standing wave light field with spatially oscillating Rabi- 
frequency. The loss rate T of a particle from the bottom of a well of the upper adiabatic 
potential U (x) into the continuum spectrum of —U(x) can be estimated semiclassically 
for a small violation of adiabaticity: 



r exp 



2? 



(Pi 



P- 



(16) 



where p± = w2m[E =p U(x)], E is the energy of the atom, and xq is determined by 
the saddle-point condition p + (x ) = p_(x ). The quantity T is the attempt frequency, 
which is of the order of the oscillation frequency u of a particle at the bottom of the 
upper potential U(x) [38j EH]- This analysis, in particular Eq. (fl6|) . can be directly 
adopted for the estimation of the Landau-Zener-type loss rate from the upper adiabatic 
potential in the simple case of one dressing field in the off-resonant scheme, and no 
dressing field in the near-resonant scheme. 

For the off-resonant scheme with only one dressing field, and for a detuning 8\ = Vq, 
the adiabatic potentials can be written as ±U(x), where U(x) = \JV(x) 2 + and 
V(x) = Vq sin(2fczr)/2. In the vicinity of the saddle point x = 0, the background 
potential can be linearised V(x) ~ V^kx, and we take the energy of the atom 
E = U(xq). We may then Taylor expand the adiabatic potentials near the saddle 
point, U(x) ~ Qi/2 + mu 2 x 2 /2, where u = 2^/V^Er/VLx with E R = h 2 k 2 /2m being the 
recoil energy in the original normal lattice. With these, one may carry out the integral 
in Eq. ffl6|) explicitly and we get for the loss rate: 



r ~ r exp 
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where the dimensionless parameter a ~ 1.32 in this case. Note that Qi characterises in 
this case the energy separation between the adiabatic potentials. Eq. (17) is different 
from the scaling relation of a typical Landau-Zener transition, where the exponent in 
the loss rate scales with the energy separation squared. One may recover this typical 
scaling relation by examining the loss rate of atoms with higher energies in the adiabatic 
potential. 

For the near-resonant lattices, different manifolds in the Floquet basis only decouple 
in the absence of dressing fields. In this case, we have only the background potentials 
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Figure 11. Landau-Zener Losses from a A/4 lattice, formed with one dressing 
transition: (a) Numerical simulation results for loss rates for atoms in the bound 
state of the higher-energy adiabatic potential in the off-resonant scheme. The solid 
line is an exponential fit to the result, (inset) The corresponding adiabatic potentials 
with parameters: Vq = 35Er, fix = 30-Er, <5i + <5 = 35Er. (b) Numerical calculations 
of the decay rate for atoms in the narrower wells in the low energy adiabatic potential 
in the near-resonant scheme, (inset) Typical potentials in the limit where the Floquet 
manifolds are weakly coupled (Vli <C f^o), with a series of Floquet manifolds each 
containing 2 dressed potentials. Parameters: f^o = 31Er, Sq = —IbEn, Six = 2Er, 
and S x = -S + (y/5% + flj + S ) x 0.8. 



±U(x) resulting from the AC Stark splitting, where U(x) = a/5q/4 + tt (x) 2 /4. We may 
also integrate Eq. ( fl~6l) to calculate the loss rate from the upper background potential: 



T ~ T exp 



—a- 



So 
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T exp 



f 3/2 
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where u' = ^VL^Er/ '5q, with a ~ 1.32. In this case, #o characterises the energy 
separation between the background potentials. 

These results suggest that the loss rate from the upper adiabatic potential can be 
made exponentially small by increasing the ratio of the energy difference of the two 
potentials and the frequency of the atomic motion in the trap. In the simple cases 
above, this translates to large Rabi-frequency Qi ^> Er in the off-resonant case, and 
large detuning S 3> E R in the near-resonant case. 

When more complicated lattice potentials are considered, the loss rate due to this 
Landau-Zener-type transition should in general be computed numerically. However, 
simple estimates of the loss rates can be obtained by assuming that for a given adiabatic 
potential the loss is dominated by the leak into the continua of the adiabatic potentials 
which are closest in energy. If we assume that the loss into the continuum of each 
potential is independent, the matrix elements for the decay into the continuum of 
any one of these nearest adiabatic potentials will have the form Eq. (|T6l) under the 
semiclassical approximation. An exponential dependence of the decay rate T on the 
energy separation between the adiabatic potentials is then obtained, similar to Eqs. (17), 
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(18). Numerically, we find the dimensionless parameter a to be on the order of unity 
or larger when the laser parameters, e.g. Rabi-frequencies and detunings, are much 
larger than the recoil energy E R . This implies that the loss rate V can always be made 
exponentially small by increasing the laser power. 

To check the general conclusion drawn above, we have also performed numerical 
simulations of the real time evolution of the atoms trapped in the adiabatic potentials by 
integrating the Schrodinger equation for the full Hamiltonian in the Floquet basis. The 
atom is initially prepared in the ground state of a potential well, and transition rates 
are extracted from the projection of the wavefunction onto the manifold corresponding 
to the initial adiabatic potential as a function of time. For the simple case in Fig. 11a, 
we vary Qi for fixed trapping frequency in the wells of the adiabatic potential ou = 
2\JVqE r /VLi. For a large gap between the potentials, the loss rate scale exponentially 
as r / ' E R oc exp(— afli/ui), as shown by the exponential fit (solid line), which is in 
agreement with the semiclassical calculation. From the fit, the factor a is estimated to 
be a ~ 0.9 and 1.2 for the exact and linearised potentials, respectively, slightly smaller 
than the results in the semiclassical approximation. In the case of near-resonant scheme, 
as in Fig. lib, we fix a/^o-^r/I^oI = 30i?R, and vary the energy gap AE between the 
adiabatic potentials. We again observe an exponential scaling, T oc exp(— a'AE/u) 
(solid line), where the trapping frequency uj ~ 15E R . In conclusion, we see one can 
always avoid large Landau-Zener-type losses from the adiabatic potential by increasing 
the laser power to make the background potential sufficiently deep. 

6. Summary and Conclusion 

To summarise, we have discussed various schemes for the generation of subwavelength 
lattices with either alkali or alkaline-earth-like atoms. The lattice spacing in these 
adiabatic potentials can be much smaller than the wavelength of the laser that creates 
the lattice potential, depending on the number of the additional dressing lasers. The 
enhanced tunneling rate one gets from the reduced lattice spacing may be useful for the 
engineering of generalised Hubbard models for quantum simulation, or for the design of 
fast quantum gates for quantum computation. The adiabatic potentials in dimensions 
higher in ID may also have interesting structures, i.e. disconnected ring potentials in 
2D, which is useful for the realisation of interesting phases. The first experimental 
steps in this direction have been recently been reported in Ref. |JT]. While we have 
discussed addressable lattices within a A/2 period, addressing on a larger scale can 
be introduced using magnetic or electric field gradients, or other recently proposed 
techniques [35], EH H21 [43] . This, combined with the gate operation, holds the prospect 
for the realisation of general quantum computation in this system. We have also 
discussed in detail the Landau-Zener-type losses from the adiabatic potentials. We 
have shown that the loss rate can always be made small at the expense of laser power. 
For large number of dressing fields, i.e. for small lattice spacing, the Landau-Zener-type 
loss would become considerable, and would be the practical limitation on how small a 



State- dependent, addressable subwavelength lattices with cold atoms 



21 



lattice spacing one can get in generating the subwavelength lattices. 
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Appendix A. Alkaline-Earth-Like Atoms and the Clock Transition 

In this appendix, we examine in detail the matrix element for the clock transition 
1 S'o — 3 -Po i n alkaline-earth-like atoms. For a detailed discussion, please refer to Ref. 

[251 ED]. 

We consider the matrix element ( 1 5'o, Imj\D\ 3 P , Fmp), where D is the electric- 
dipole operator. In the presence of a hyperfme interaction, \ 3 Pq, Ftuf) is not an 
eigenstate of the system. Under the first order perturbation, the correction to the 
hyperfme sublevels in 3 P manifold can be written as: 

v^i t 17 ,{l,J,Fm F \Vu\ z Pa,Fm F ) 

x> j ^> — m - £(3p o) — > (A - 1} 

where 7 labels the manifolds that the the state is coupled to, and Vm is the hyperfme 
interaction. Due to the scalar nature of Vhf, the total angular momentum F and its 
projection on z-axis nip are good quantum numbers, as is shown in the expression 
above. As we are only interested in the El transitions from x Sq to 3 Po, we only consider 
the nuclear magnetic dipole contribution in the hyperfme interaction, which fixes J = 1 
in the expression above. 

The matrix element for the clock states then becomes: 

CSaJmilDl^Fmp) 

= ^fSo, Im I \D\'y, J —I, Fm F )A 

7 

= E C^^SoImjlDl^mm^A, (A.2) 

mj+mj=mF 

where A = ^'^g^gpp^"*^ , and C^ ljmi is the Clebsch-Gordon coefficient for the 
state I7, J = 1, F). 

For the convenience of discussion, we now consider the concrete example of m Yb, 
which has a nuclear spin of I = 1/2. It is obvious that the matrix element is nonzero for 
both the 7r-polarised laser field and the circularly-polarised laser field. For 7r-polarised 
lasers, the matrix element is effectively given by 

(7, J = l,mj = 0,m 7 = ±-\D\ 1 S ,m J = 0,m 7 = ±-); (A. 3) 

while for circularly-polarized lasers, the matrix element is effectively given by 

(7, J = l,mj = ±l,m 7 = T^\D\ l S ,m,j = 0,mj = ( A - 4 ) 
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for the cr + - and a~- polarised light, respectively. 

The ratio between the matrix elements in the two cases is given by that of the 
corresponding Clebsch- Gordon coefficients in Eq. (A2). One may further relate these 
matrix element to the natural linewidth T of the clock transition: 

^0 

where Uq is the frequency for the clock transition. The constant factor (3 comes from the 
Clebsch- Gordon coefficient, which has j3 = 1/3 for 7r-polarised lasers and (3 = 2/3 for 
circularly polarised lasers. For 171 Yb, r ~ 2n x lOmHz [271 f28l [29] . and a typical laser 
intensity of I = 10kW/cm 2 gives a Rabi-frequency of ~ 2tt x 220kHz for a circularly 
polarised field. 
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